function c = fn_c(y,theta,T)

c = y + theta.*T;

end